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Preliminary  Note 

Application  of  the  Mean  Spherical  Approximation  to  the  Estimation  of  Single  Ion 
Thermodynamic  Quantities  of  Solvation  for  Monoatomic  Monovalent  Ions  in 

Aqueous  Solutions 

W.  Ronald  Fawcett 
Department  of  Chemistty 
University  of  California 
Davis,  California  95616,  U.S.A. 
ai*4 

Lesser  Blum 

Department  of  Hiysics,  P.O.  Box  23343 
University  of  Puerto  Rico 
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Introduction 

The  separation  of  thermodynamic  quantities  ftv  electrolytes  in  sdution  into  contributions  fix* 
the  cation  and  anion  is  a  well  known  problem  in  solution  electrochemistry  [1].  It  has  inqxmant 
consequences  widt  respect  to  the  establishment  of  an  absolute  potential  scale,  a  sul:gect  whidi  has 
been  considered  in  detail  by  Trasatti  [2, 3]  and  Parsons  [4].  Obviously,  an  extiadietmodynainic 
assunqrtion  is  required  to  achieve  this  separation  and  various  approaches  based  on  approptiiue 
theories  have  been  considered  [1-4]. 

Way  recently,  we  discussed  the  applicadmi  of  the  mean  spherical  appioximatioo  (MSA)  to 
estimate  die  (Hbbs  energies  of  solvation  of  monoatomic  monovalent  ions  in  ptdar  atdvents  [5]. 
According  to  the  MSA  model  two  solvoit  dependent  i»tameten  other  thra  those  used  in  die 


classical  Bom  model  are  required.  These  are  the  polarization  parameter  X  which  accounts 
especially  for  shOTt  range  intoactions  between  the  solute  ion  and  surrounding  solvent  dipoles,  and 
rs  the  radius  of  the  solvent  molecule  represented  as  a  sphere.  It  was  shown  that  the  ratio  XAs 
depends  on  whether  the  tmnioatomic  ion  is  an  alkali  metal  ion  or  a  halide  ion,  and  cm  the  nature  of 
the  solvent  such  that  it  increases  with  solvent  basicity  in  the  case  cations,  and  widi  solvent 
acidity  in  the  case  of  anions  [5]. 

The  purpose  of  the  present  paper  is  to  examine  father  the  abiliQ'  of  the  MSA  model  to 
estimate  single  ion  thermodynamic  quantities  and  thus  evaluate  its  significaiK:e  as  an 
extrathermodynamic  route  to  these  quantities  for  simple  monoatomic  ions.  In  tins  regard,  it  is 
important  to  have  reliable  estimates  of  the  radii  of  the  ions.  These  parameters  are  most  often  based 
on  the  values  given  by  Pauling  [6]  which  were  estimated  using  crystallognqthic  data.  However, 
in  recent  years,  extensive  information  has  bectxne  available  regarding  the  structure  of  aqueous 
electrolyte  solutions  on  the  basis  of  neutron  and  X-ray  difftaction  experiments.  This  infionaatkm 
has  been  compiled  by  Marcus  [7]  and  provides  an  addititmal  source  of  ionic  aze  data  wfaidi  can  be 
compared  with  those  proposed  by  Pauling  [6].  The  ttifects  of  tiie  ionic  radii  mi  die  estimation  of 
Gibbs  solvation  energies  for  single  ions,  and  the  resulting  MSA  parameters  are  described  in  die 
present  note. 

The  Model  and  Its  Application  to  Thermodynamic  Data 

On  the  basis  of  the  MSA  [S,  8-10]  the  standard  Gibbs  energy  of  ion,  i,  in  solvent,  s,  at 
infinite  dilution  is  given  by 


qO  (jj  _  .  / j .  JL'l _ 1 

V  ®sj 


8n  eo 


(1) 


(n+Sj) 

where  No  is  the  Avogadro  constant,  eo,  the  electronic  charge,  zi,  the  ionic  valence,  Cq,  die  _ 

r 

permittivity  of  free  space,  £$,  the  dielectric  constant  of  the  pure  solvent,  n,  die  ionic  radius,  and  5s,~ 
the  MSA  parameter  which  depends  on  the  nature  of  the  solvent  and  whether  die  km  is  a  cation  or 

anion.  The  last  quantity  is  related  to  the  polarization  parameter,  X,  and  stdvent  radius,  rg, - 

discussed  above,  by  die  equation 


Avatlatlllty  Codas 
Avail  and/or 
Dial  Spaoial 


□  □ 


8s  =  Tj/X  (2) 

In  earlier  versions  of  the  MSA,  the  quantity,  X,  was  estimated  from  the  Werdieim  equation  [11] 
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that  is,  using  the  relationship 

X2(l+X)^  =  16es  (3) 

However,  the  resulting  estimates  of  X  do  not  account  for  inqtortant  parts  of  die  short  range 
electrostatic  interactions  between  die  ion  and  surrounding  dipoles,  that  is,  die  "stickiness"  of  ion- 
dipole  interactions.  These  interactions  result  from  higher  wdcr  moments  associated  with  the 
electrostatic  description  of  the  solvent  and  from  its  polarizability.  In  our  recent  analysis  of  Gibbs 
solvation  data  fw  monovalent  ions  [5],  X,  and  thus  8s,  were  treated  as  adjustable  parameters 
determined  from  single  ion  data  reprated  in  die  literature  [12, 13].  In  die  present  paper  we  af^ly 
this  analysis  to  thermodynamic  data  for  the  alkali  metal  halide  salts  in  water,  ft  should  be  noted 
that  the  quantity  estimated  by  eq.  (1)  is,  stricdy  speaking,  only  the  uxi-d^le  portion  of  the  Gibbs 
solvation  energy.  Thus,  it  is  assumed  in  the  present  treatment  that  one  can  avoid  eiqididdy  treating 
contributions  from  dipole-dipole  and  shrat  range  electrostatic  interacticHis  by  regarding  5b  as  an 


adjustalde  parameter  whose  value  accounts  for  these  effects. 

The  exinession  for  the  standard  Gibbs  enogy  of  solvation  of  die  mmovalent  salt  CA 
becomes 


g;(CA) 


> .  jw 

Sk  Co 


(4)[: 


rc+8sc 


(4) 


where  rc  and  ta.  are  the  radii  of  the  cation  and  anion,  respectively,  and  8sc  and  8ba>  die 
ctffiesponding  values  of  the  MSA  distance,  8$.  It  is  ^parent  that,  if  one  knows  the  values  of  the 
ionic  radii,  only  two  adjustable  parameters  namely,  8sC  and  5sa>  are  involved  in  estimating  die 
Gibbs  solvation  energy  for  the  salt  This  expression  was  used  by  Latimer  et  aL  [14]  to  account  for 
the  variation  in  the  Gibbs  energy  of  solvatitxi  of  the  alkali  metal  halides  with  size  of  the  congiooent 
ions  but  the  interpretation  of  die  parameters  8sc  and  8sa  was  quite  different  They  were  regarded 
as  corrections  to  the  Pauling  radii  necessary  so  that  the  Bmn  model  could  be  ^iplied  to  the 
problem.  Anotherwayrtf  interpreting  these  quantities  with  the  context  of  a  primitive  model  is  as 


effective  distances  over  which  the  dielectric  constant  varies  frcnn  a  low  value  near  the  ion  to  die 
bulk  value  furtho’  away  [IS]. 
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The  entn^y  of  a  single  ion  may  be  found  firom  the  temperature  derivative  of  the  Gibbs 


eneigy: 


®  Sjcco  W  ^  (n+8s)  I  esj  (r,+( 


1  d^‘ 

(n+6s)2 


Its  estimation  requires  two  additional  parameters  namely  the  temperature  coeffident  of  the  didectric 
constant  de^dT  which  is  available  from  experiment,  and  the  teii^)erature  coefficient  of  die  MSA 
parameter,  dVdT.  The  first  term  on  the  right  hand  side  of  eq.  (S)  may  be  regarded  as  giving  the 
contributitm  to  the  entr(q}y  due  to  long  range  electrostatic  interactions,  and  die  seccmd,  diat  due  to 
short  range  interacticms.  In  earlier  wmk  [14],  d»  parameter  8s  was  ctHisideted  to  be  a  constant  so 
that  an  inqiortant  contribution  to  the  entropy  was  ignmed.  >^thin  the  context  of  the  MSA,  the 
temperature  variation  in  5s  is  due  to  a  corresponding  change  in  dVdT  (eq.  (3)).  Since  the 
"stickiness"  of  ion-solvent  interactions  are  expected  to  decrease  as  tenqierature  increases,  dX/dT  is 
expected  to  be  a  negative  quantity  and  dS^dT,  positive.  For  the  same  reas(»s  dut  die  relationship 
between  X  and  es  is  considered  to  be  more  complex  than  diat  given  by  eq.  (3),  the  relationship  that 
could  be  derived^fitom  this  equation  between  dX/dT  and  dEs/dT  does  not  properly  account  for  the 
change  in  ion-solvent  interactions  with  temperature. 

It  is  convenient  to  re-express  the  single  ion  entropy  in  terms  of  the  single  km  Gibbs  energy. 
Thus,  eq.  (5)  may  be  rewritten  as 


^s^*^  dfis  8x  Co 


/'i  ^  >  s' '  dCs  8x  Co  ._o  r.vo 


When  estimates  of  G°  (i)  are  available,  entrc^ic  data  for  a  series  of  ions  may  be  {dotted  on  die 
basis  of  eq.  (6)  to  obtain  a  value  of  the  parame^  dSsAlT. 

Experimental  values  of  the  Gibbs  energies  of  solvation  for  20  alkali  metal  halide  aalta  based 
on  the  five  cations  finrn  Li'*’  to  Cst  and  die  four  anions  from  P  to  I*  were  obtained  by  aAHng  die 
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corresponding  single  ion  values  from  the  tabulation  by  Parsons  [4].  The  single  ion  values  are 
those  based  on  the  assumption  that  (H'*')  =  0.  and  should  give  the  best  estimates  of  the  required 
electrolyte  data.  The  values  of  (CA)  for  the  20  salts  were  then  fitted  to  eq.  (4)  using  non-linear 
least  squares  [16]  to  extract  the  best  values  of  5sc  and  Sja  for  a  given  set  of  ionic  radii.  Quite  good 
fits  of  the  experimental  data  to  the  MSA  equation  for  the  Gibbs  energy  solvatkm  of  the  1-1 
electrolytes  formed  fiom  these  simple  itrns  were  obtained  for  both  sets  of  radii  However,  cm  die 
basis  of  the  estimates  of  the  Gibbs  solvation  energies,  the  radii  obtained  fiom  solutitm  data  are 
cleariy  preferred  since  the  values  of  G°  (i)  are  very  close  to  experimental  values  of  the  real  solvation 
energy  for  single  ions  [17,18],  also  tabulated  by  Parsons  [4].  Of  course,  this  quanti^,  which  is 
equal  to  the  negative  value  of  the  work  function  for  the  ion  [19],  contains  a  contributkxi  fiom  die 
surface  potential  of  water  which  is  estimated  to  be  approximately  130  mV  at  25  X  [4, 20].  As  a 
result,  the  experimental  real  potential  is  about  12.5  kJ  mol'^  less  negative  than  die  single  Gibbs 
energy  of  solvation  for  cations,  and  about  12.5  kJ  rool*^  more  negative  for  anions.  Considering 
the  analytical  simplicity  of  the  MSA  expression  for  the  Gibbs  solvation  energy,  the  present  results 
are  indeed  very  gratifying.  Using  the  single  ion  data,  (me  may  estimate  a  value  for  the  constant 
required  to  ccmvert  the  data  on  the  scale  referenced  to  G^  (H'*') »  0  to  an  absolute  scale.  In  the  case 
of  the  radii  based  on  diffracticm  data,  the  resulting  constant  is  -1091  ±  7  kJ  mol'^  which  ccmpues 
well  with  the  best  value  for  the  real  potential  of  the  H+  ion  (-1088  ±  2  kJ)  [17, 18].  Thus,  the 
results  reported  here  provide  an  acceptable  basis  for  calculating  single  electrode  potentials. 

The  present  analysis  was  carried  further  by  examining  the  data  for  the  entrt^  of  sotvaticm  (tf 
the  same  ions  in  aqueous  solutions  using  eq.  (6).  Tabulated  values  fyc  this  quantity  are  also 
referenced  to  the  solvation  entropy  of  the  hydrogen  ion  (S®  (H+)  m  0)  [4].  Keeping  diis  in  mind, 
eq.  (6)  can  be  rewritten  as 
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Where  Y  =  AS;(i)  (l  -  M  (8) 

V  ®s/  ®s 

and  AS°  (i)  is  the  tabulated  entropy  of  solvation  referenced  to  that  for  km.  The  variables  X 

and  Y  were  calculated  using  the  best  values  of  the  Gibbs  solvation  enerf^es  firom  Table  1.  A  [doc 
of  the  data  for  the  alkali  metal  cations  is  shown  in  Fig.  1,  an  excellent  linear  oondatioa  being 
obtained  (r  »  0.993).  The  resulting  value  of  (H'**)  estimated  firom  die  inteicqpt  is  -126 12  J  K*^ 
this  agrees  quite  well  with  die  value  estimated  by  Conway  (-131 J  K*i  mol*^)  [1].  The 
value  of  dSgc  /  dT  obtained  from  die  slope  is  0.029  pm  K*^  A  more  important  test  of  die  iffesoit 
model  is  the  comparison  with  results  obtained  firom  the  analysis  of  data  for  amoas.  The 
corresponding  plot  for  the  halide  ions  is  also  ^own  in  Fig.  1.  The  linear  oonelatioo  in  dds  case  is 
not  as  good  (r  =  0.988),  and  the  estunace  of  is  - 1 12 1 3  J  K**  nxd*^.  Keqnng  in  mind 

that  data  are  available  for  tmly  four  ions  this  result  agrees  reasonably  well  widi  dist  based  on  die 
alkali  metal  halide  data.  Estimates  S°(H^)  based  on  die  Pauling  radii  are  tether  finm  the 
estimate  of  Conway  [1]  and  farther  from  each  other  when  one  compares  cadooic  and  anionic  dt»fa; 
therefore  they  are  considered  less  reliable.  Hnally,  the  value  of  dfi^A/  dT  from  the  skpe 

of  the  plot  for  the  anionic  data  is  0.047  pm  K*^.  Just  as  the  values  of  5a  are  difteent  fior  cations 
and  anions,  so  are  dieir  temperature  coefficients. 

Discussion 

The  above  analysis  demonstrates  that  the  MSA  provides  an  excellent  desciqition  of  die 
thermodynamic  pnperties  for  simple  monovalent  ions  in  aqueous  solution.  Having  dioaen  a  set  of 
ionic  radii,  only  two  additional  adjustable  parameters  are  required  to  achieve  a  fit  with  experimemal 
data  for  the  Gibbs  solvation  energy,  namely,  the  MSA  distances  fix  cations  and  •nirni*.  It  should 
be  emphasized  duu  die  values  of  diis  parameter  repxted  here  are  quite  difiteent  finxn  diat  obtained 
fincrni  the  Wertheim  equatkm  which  gives  an  estimate  of  S1.6  pm  fix  water  at  25  ^  when  the 
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molecular  diameter  is  assumed  to  be  274  pm.  As  one  would  expect,  small  diffoences  in  8$  lead  to 
lai^e  changes  in  the  estimates  of  (i).  These  differences  are  attributed  to  die  fact  that  the  original 
Wertheim  model  represents  the  pure  solvent  as  a  collection  of  hard  spheres  with  point  dipoles. 
Although  estimation  of  5$  on  the  basis  of  eq.  (3)  accounts  for  some  of  the  other  electrostatic 
interactions  between  solvent  molecules,  this  value  does  not  account  for  the  mote  con^lex  aspects  oi 
ion-dipole  interactions  in  electrolyte  solutions.  This  is  most  easily  seen  when  one  recognizes  diat 
the  solvation  of  cations  and  anions  with  equal  radii  is  not  the  same.  Thus,  the  present  empirical 
approach  to  dealing  with  the  "stiddness"  of  these  interactitxis  seems  appropriate  when  ont 
considers  the  complexity  of  the  theory  that  would  be  required  to  provide  an  adequate  descr^tion 
including  higher  order  electrical  moments  and  polarizalnlity. 

It  is  also  apparent  from  the  present  analysis  that  die  choice  of  ionic  radii  used  in  estimating  die 
thermodynamic  properties  is  extremely  important  It  is  also  not  surprising  that  the  best  fit  between 
theory  and  experiment  was  obtained  using  the  radii  from  diffraction  experinoents.  In  our  previous 
analysis  of  the  data  fcx*  the  Gibbs  solvation  energy,  the  single  ion  values  estimated  by  Cc»way  [1] 
and  Pauling  radii  were  used  to  obtain  the  best  fit  between  the  experimental  data  and  the  MSA  model. 
However,  the  present  method  of  separating  the  cationic  and  anionic  components  using 
thermodynantic  data  for  the  salts  provides  an  extradtermodynamic  route  to  the  sin^  ion  data  and 
produces  results  which  are  much  closer  in  agreement  with  the  real  solvation  energy  data  of  Randles 
[17].  The  uncertainty  in  most  radii  used  to  obtain  this  separation  is  in  the  order  of  ±10  pm.  As 
more  data  are  gathered  from  the  dif&action  experiments,  the  resulting  estimates  of  q  should  inqnove 
and  may  provide  a  better  description  of  the  single  ion  quantities. 

In  order  to  estimate  the  entropy,  and  eventually  the  enthalpy  of  solvation,  two  more 
parameters  are  required,  namely,  the  temperature  coefficients  of  the  MSA  distances  for  cations  and 
anions.  The  present  analysis  requires  that  one  uses  estimates  of  the  single  ion  Gibbs  solvatkm 
energies  to  obtain  these  parameters  and  yields  excellent  results.  The  major  point  of  weakness  m  die 
analy^  is  uncertainty  in  the  values  of  q.  Thus,  inqirovement  in  these  estimates  should  lead  to 
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results  from  the  entropic  analysis  which  agree  widi  one  another  when  one  cmiqMires  data  for  cations 
and  anions. 

The  pr^nt  analysis  for  tlw  aqueous  solvation  data  can  easily  be  ectended  to  other  pdar 
solvents  using  data  fear  the  Gibbs  energy  and  entropy  of  transfer  of  simple  1-1  dectrolytes 
composed  of  aiifaH  metal  cations  and  halide  anions.  Unfortunarely,  diese  data  are  usually  available 
for  less  than  ten  electrolytes  from  the  group  considered  here  [12, 21J.  This  means  dtat  die  chances 
of  obtaining  an  unambiguous  fit  using  non-linear  least  squares  are  smaller.  However,  such  an 
analysis  would  be  very  helpful  mainly  because  it  would  allow  one  to  assess  die  TATB  assumption 
which  has  been  piedominandy  used  to  separate  cationic  and  anionic  conqionrats  for  themxxfynamic 
transfer  data  [12, 13].  This  assumption  is  based  on  equal  solvation  of  large  ions,  namdy,  die 
tetrqihenylarsonium  cation  and  the  tetraphenylberate  anion.  On  the  basis  of  die  MSA,  etpial 
solvatimi  of  these  ions  would  require  not  only  that  their  effective  radii  be  equal,  but  also  diat  dieir 
MSA  parameters  be  the  same.  Since  these  ions  are  tetrahedral  and  not  q[)heiical,  the  stated 
conditions  may  be  met  Separation  of  the  cationic  and  anionic  conqxments  of  thermodynamic 
parameters  for  simple  electrolytes  on  the  basis  of  die  present  model  would  provide  an  aocqitable 
alternative  by  which  the  TATB  assumption  could  be  tested.  Tliis  subject  will  be  dealt  widi  in  detail 
in  a  future  paper. 
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Table  1.  Values  of  Ionic  Radii,  MSA  Parameter  Ss  and  Estimated  Gibbs  Energy  of  Sdvadon  for 
the  Alkali  Metal  Cations  and  Halide  Anions 

Ion  Pauling  Dilfiacdon  Data 


r,/pm 

-Gj(i)/kJmoH 

r,/pm 

-Gj(i)/kJmoH 

u+ 

60 

481 

71 

503 

Na+ 

95 

386 

97 

423 

K+ 

133 

318 

141 

332 

Rb+ 

148 

298 

150 

319 

Cs+ 

169 

273 

173 

288 

F- 

136 

445 

124 

432 

a- 

181 

345 

180 

320 

Br 

195 

322 

198 

295 

I- 

216 

293 

225 

264 

MS  A  Parameter 

5s/)pm 

Cadcms 

82.4 

65.3 

Anions 

18.0 

34.6 

Standard  Deviadon 
for  fit  to  eq.  (4) 

1.5  U  mol-i 

l.lkJmol-i 

J 


Legends  for  Figures 

Figure  1.  Plot  of  entropy  of  solvation  corrected  for  long  range  electrostatic  effects,  Y  (eq.  8) 
against  X,  a  quantity  proportional  to  the  square  of  the  Gibbs  solvation  energy  (eq. 
9)  for  the  alkali  metal  and  halide  ions  in  aqueous  solutions. 


